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Abstract. A new "Dynamical Mean-field theory" based approach for the Kondo 
lattice model with quantum spins is introduced. The inspection of exactly 
solvable limiting cases and several known approximation methods, namely the 
second-order perturbation theory, the self- consistent CPA and finally a moment- 
conserving decoupling of the equations of motion help in evaluating the new 
approach. This comprehensive investigation gives some certainty to our results: 
Whereas our method is somewhat limited in the investigation of the J < 0- 
model, the results for J > reveal important aspects of the physics of the 
model: The energetically lowest states are not completely spin-polarized. A band 
splitting, which occurs already for relatively low interaction strengths, can be 
related to distinct elementary excitations, namely magnon emission (absorption) 
and the formation of magnetic polarons. We demonstrate the properties of the 
ferromagnetic Kondo lattice model in terms of spectral densities and quasiparticle 
densities of states. 



PACS numbers: 71.10.Fd, 75.30.MB,75.30Vn 

1. Introduction 

The Kondo model and its periodic extension, the Kondo lattice model (KLM), which 
describe spin-exchange interaction between a localized spin or a system of localized 
spins, respectively, and a band of itinerant electrons, has been subject of intense 
theoretical studies in the past [Q, ||, ||, ^, ^, ^. This model has been applied to a 
variety of different problems in solid-state physics using both a ferromagnetic and 
antiferromagnetic coupling constant J. 

The model with J < is the one originally known as Kondo lattice model or simply 
Kondo model in its non-periodic form with a single impurity spin in the system. It 
was used by Kondo to explain the unusual temperature behavior of the resistivity of 
magnetic impurities in non-magnetic hosts Q . The negative spin-exchange interaction 
can be derived from the hybridization of a correlated "atomic" level with a conduction 
band, the situation described by the Anderson model iQ, ||. In the limit of a low- lying 
half- filled atomic level and strong correlations, the Anderson model can be mapped 
onto the Kondo model with a negative exchange constant ^ . The Kondo lattice model 
is still subject to much theoretical work, the main objective is the understanding of 
the unusual physical behavior found in Heavy-Fermion materials jsj. 
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A model with identical operator structure in the Hamiltonian, but with positive 
exchange constant has been known in the literature for a long time by many different 
names (double exchange model, s-d model, s-f model,...) |jl|, ^, ||, ^. For clarity, 
we will refer to this model in the following as ferromagnetic Kondo lattice model. 
The model with ferromagnetic exchange has to be understood as an effective one. 
The origins of the exchange with J > are found in the interband Coulomb 
correlations This situation is believed to dominate the physical properties of 
important systems such as the magnetic semiconductors ||l^ (EuX; X = O, S, Se, 
Te), the diluted magnetic semiconductors (Cdi_a;Mn^Te, Hgi-^^Fea^Se), and the 
"local moment" metals (Gd, Dy, TbV To these problems, the ferromagnetic 
KLM was successfully applied ^ Recently, this variant of the KLM 

has gained a lot of interest with the discovery of the colossal magnetoresistance 



(CMR) materials |16, |17|. In these materials, typically manganese oxides with 
perovskite structure (Lai_a;(Ca,Sr)a;Mn03), the double-exchange model |^ has 
been successfully applied to explain the origin of ferromagnetic order and is expected 
to be a good starting point to investigate the resistivity anomalies This 
double-exchange model, however, is nothing else than the Kondo lattice model with 
ferromagnetic (positive) exchange constant in the strong coupling limit. In the 
CMR materials, the localized S = |-spin of the model represents the more or less 
localized manganese 3d-t2g electrons, whereas the conduction band is formed by the 
Cg electrons. The interband-exchange interaction is nothing else but the intra-shell 
Hund's rule coupling. Since the Sd-e^ electrons of the manganese form a relatively 
narrow band (theoretical results from band-structure calculations: 1 — 2eV |2^, pH 
and experimental estimates: 3 — AeV p2| , |2^ ) and Hund's coupling is assumed to be 
large, the model has to be taken in the intermediate to strong coupling regime. There 
are few estimates about the value of the interaction constant in the literature, e.g. 
J leV ||l^, but these are challenged as to be too small Most theoretical 
papers of the last years concerned with colossal magnetoresistance assume classical 
spins 5 — > oo ||2^, 25, This has been justified by the assumption of 



JS oo |25|. Although it is true that the important energy scale is JS, there are 
much more implications of — > oo that are not justified in the strong-coupling limit 
for a S" = I system. In several papers, it was stated that ". . . the Cg electrons are 
oriented parallel to the t2g spins." or equivalently ". . . so one only need consider 
configurations with Cg electrons parallel to core spins." H^. We will show below using 
exact results as well as several well-defined approximation methods, that for S = ^ 
there is a considerable amount of spin-| spectral weight located in the main region of 
the spin-t states even for large interaction strengths. The assumption of a half-metallic 
state [Q, made in the two citations above can therefore never be met in the KLM 
with quantum spins and is merely an effect of the (unphysical) limit of "classical" 
spins. The recently discussed half-metallic behaviour of the manganites |^ must 
have a different origin. 

However, for the opposite sign of J, exactly the assumed effect happens in the 
strong-coupling limit: the lowest-lying excitations in the conduction band density of 
states will be purely spin-|. This already implies that results for the Kondo lattice 
model with J > and J < cannot simply be reverted into the respective other case. 
The change of sign changes the whole physics of the system. For J < an antiparallel 
( "antiferromagnetic" ) alignment of the conduction band spin and the localized spin 
lowers the internal energy. For a sufficient band filling, this tends to a screening of 
the local moments by conduction electrons, well-known from the Kondo effect that 
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refers to a single magnetic impurity in a conduction electron sea. From this, the name 
"Kondo lattice model" was originally derived for the J < case. 

We will further show that already for comparatively low interaction strengths the 
spin-exchange interaction alone leads to an opening of a gap in the density of states. 
This extraordinary correlation effect could give hints to the explanation of a recently 
discovered pseudogap in the managese oxides | |3^ . 

To prove our claims already laid out so far, we will first review two important non- 
trivial exactly solvable limiting cases of the Kondo lattice model in Sec. ^ The first 
is the zero-bandwidth limit ("atomic limit") where the bandwidth of the conduction 
band is set to = js^]. The second exactly solvable limiting case is the so-called 
ferromagnetically saturated semiconductor |3^, |3^, 39 , ^ ^ ^ . 



This is essentially the zero-temperature limit of the model with vanishing electron 
density and fully aligned spin system. In this limit, striking correlation effects can 
be observed and discussed. These limiting cases will already give clear evidence to 
our propositions made above. In Sec. ^ we will present a new dynamical mean-field 
theory (DMFT)-based approach for the KLM with -S* = ^ spins. To circumstantiate 
our theory, we will introduce three more approximation schemes which also keep the 
spin as a quantum variable, not relying on the classical spin limit. The first will be a 
second-order perturbation theory (SOFT) for the KLM based on the projector operator 
formalism |43, Q. The self- consistent CPA (coherent potential approximation) is a 



straightforward extension to the well-known CPA for the KLM |, |6[ 
It starts from the zero-bandwidth limit discussed before. The third approximation 
method is a moment-conserving decoupling procedure for the equation of motion of the 
single-electron Green function. This approximation scheme continuously evolves from 
the exactly solvable limit of the ferromagnetically saturated semiconductor. It will 
be called Moment-Conserving Decoupling Approximation (MCDA). The comparison 
of the results obtained by these three methods and our DMFT scheme, each of which 
starts from a different limit allows to evaluate the range of applicability of the new 
approach, and to select the most trustworthy common features of all methods to gain 
a reliable picture of the physics of the ferromagnetic KLM. 

2. The Kondo lattice model and its many-body problem 
2.1. Hamiltonian 

The Kondo-lattice model (or s-f model) traces back the characteristic features of the 
underlying physical system to an interband exchange coupling of itinerant conduction 
electrons to (quasi-) localized magnetic moments described by the following model 
Hamiltonian ^ 

H = H,+H,f + iHu + Hff)= (1) 



^irj Sfc '^Io-^~*'°^') ^"^^ ^i<^ respectively, creation and annihilation 

operators of a band electron being specified by the lower indices. The hopping 
integrals Tij are connected by Fourier transformation to the single-electron Bloch 
energies r„ = ^ Efe e(fc)e*'^ (^*--^^). 
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The interband (s/) exchange with couphng strength J is taken as an intra-atomic 
interaction between the conduction electron spin ai and the locahzed magnetic moment 
represented by the spin operator Si. For practical reasons it is sometimes convenient 
to use the second quantization representation of the band electron spin ai which leads 
to the following form of the interband-interaction term: 



^ - V E f^-^^"- + S^cl^C^a) (2) 



Here we have used the abbreviations Uicr = c\^Cia, = -l-l(— 1) and S'f = 

S'f + iz„Sf . The first term in (||) describes an Ising-like interaction between the 
z-components of the spin-operators, while the second term incorporates spin exchange 
processes between the localized and the itinerant system. 

The last two terms are an extension to the original model: Hu leads to the 
'correlated Kondo lattice model', it introduces correlations between the conduction 
electrons in form of a 'Hubbard-type' interaction. We will include this term in 
some parts of the discussion below. The second term, Hff represents a direct 
spin-exchange between localized moments on different lattice sites. Although the 
sf exchange can lead to an effective RKKY interaction between the spins for non- 
vanishing band occupation, a "superexchange" could also be included. In case of an 
empty conduction band, this term becomes essential as source of magnetic order. Jij 
are the superexchange integrals. We state once more that the original Kondo-lattice 
model (or s-f model) is defined by H = Hg + Hgf, only. The additional terms in 
Eq. (|l|) are used as soon as physical requirements do not allow to neglect them. 

If we are mainly interested in the conduction electron properties then the single- 
electron Green function Gija{E) — {{cia-', Cj^)) e is of primary interest. Its equation of 
motion reads for the correlated KLM 

{E6i,n — Tijn)Gmja{E) = HSij — - J {ZcrIiiJa{E) + Fn^j„{E)) + UTiiijcr{E) (3) 



E 



The two types of interaction terms in let appear the "spinflip function" Fim.ja (E) = 
{{S~'^Crn^a-,c^j„))E and the "Ising function" Iim,3a{E) = {{SiCm.a]c]^))E, while the 

"Hubbard-function" Tihn.ja{E) = {{c\_^ci^„Cma;c^jcr)) e comes into play only when 
the "Hubbard-interaction" (last term in Eq. (|l|)) is switched on. 

The three "higher" Green functions on the right-hand side of Eq. (||) prevent 
a direct solution of the equation of motion. A formal solution for the Fourier- 
transformed single-electron Green function, 

GUE) = ((c..;cL)). = E_,^^)\^^^^E) (4) 

defines the electronic selfenergy J^k^iE) by the ansatz {{[cka', El — Hs]_ ; c]^^))e = 
'^ka{E)Gka{E) For the general case Y,kcr{E) cannot be determined rigorously. 

Before introducing some approaches to the not exactly solvable many-body 
problem of the KLM let us discuss in the next sections two rather illustrative limiting 
cases which can help to test the unavoidable approximations. 

2.2. The zero-bandwidth limit 

Let us assume that the arbitrarily filled conduction band is shrinked to an A^-fold 
degenerate level Tq: e(fe) Tq Vfc. Nevertheless, we consider the /-spin system as 
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collectively ordered for T < by any direct or indirect exchange interaction. In this 
case, the hierarchy of equations of motion decouples exactly and can rigorously be 
solved . The resulting energies and respective spectral weights are 

Ei=To- ^JS ai, = ^^^{S + 1 + z,{S') + 1-a - ^.A-. - {S + 

^2 = To + i J(5 + 1) a2a = {S - z„{S') - 7_, + z,A_, - S{n^^)} 

Es^n + U- i J(5 + 1) - {S{n^,) - 7-. + Za^-a} 

Ei^To + U+ ^JS - {{S + + - ^,A_,} 

The "Hubbard-U" in and E4 indicates that these excitations are bound to a 
double occupancy of the respective lattice site. Ei and E2 appear when our "test 
electron" enters an empty site. If it orients its spin parallel to the local /-spin then 
the energy Ei is needed. In case of an antiparallel spin orientation a triplet or a 
singlet state is formed. The first requires the energy Ei, the second £^2- The latter 
is therefore two-fold degenerate. The spectral weights are, contrary to the energy 
levels, strongly dependent on the magnetization state of the / system and the band 
filling. For a complete solution one needs the average occupation number {n_„) and 
the mixed correlation functions 70- = {Sf c^^_^Cia) and A_o. — {Sfnia). The evaluation 
can selfconsistently be done by use of the spectral theorem for the Green functions 



Giia{E), Iii^ia{E) and Fn^i„{E) (cf. Ref. |33 ). It is interesting to observe that in any 



case from the four poles only three do appear. For less than half-filled bands and 
J > (J < 0) aso- {uia) vanishes, and for more than half-filled a2a (aio-) does. It 
should be mentioned that the spectral weights aia do not explicitly depend on the 
coupling constants J and U . That means, on the one side, that even for U — Q the s-f 
interaction produces a splitting into four not coinciding quasiparticle levels. On the 
other hand, there is a striking dependence of A_cr and {n-„) on the sign of the s-f 
coupling. That transfers to the spectral weights giving them some indirect dependence 
on J/\J\- For J > and J < the order of the energy levels is different resulting 
via the spectral theorem in different correlation functions. Note that the mentioned 
dependence on J concerns only the sign of J . The spectral weights are not influenced 
by the absolute value |J|. 

2.3. The ferromagnetically saturated semiconductor 

There is another very instructive limiting case that can be treated rigorously. It 
concerns the situation of a single electron in an otherwise empty conduction band at 
T — 0, when the local moment system is ferromagnetically saturated. In this case the 
Coulomb interaction is meaningless, the "Hubbard-function" TumjaiE) is identical 
to zero. In the zero-bandwidth limit, discussed in the last section, for the |-spectrum 
all spectral weights disappear except for ai] = 1, while for the |-spectrum the levels 

El and E2 survive with aix = jE+i ^^'^ '^'21 — 2S+1 • 

For finite bandwidth the mentioned special case is that of a ferromagnetically 
saturated semiconductor ||, |3|, || ||, ||, |3|, |o[ • 1^ this case, the 

spin-f quasiparticle density of states p^{E) is only rigidly shifted compared to the 
"free" Bloch density of states. 

p^^E)^-^^%^^E+\jS) (6) 
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Figure 1. Spin-J. spectral densities (quasiparticle bandstructures) of the 
ferromagnctically saturated semiconductor for •S' = ^ and J as indicated along 
several symmetry directions. The left (right) column displays the J < (J > 0) 
case, the top-most picture shows the interaction-free spectral density. 



Consequently, the quasiparticle dispersion is undeformed with respect to the Bloch 
energies, E^{k) e{k) — \ JS. The J,-spectrum is more complicated because a J,- 
electron has several possibilities to exchange its spin with the antiparallcl / spins. 
Therefore, the a =J,-spinflip function is not at all trivial. However, its equation of 
motion decouples exactly, producing a closed system of equations which can be solved 
after Fourier transformation for the single-electron Green function. The corresponding 
selfenergy llki{E) reads: 

^''(^) =^E(^-^(^-9) + ^-^^-M9)) (8) 
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hw{q) are the spin wave energies following from the Heisenberg exchange Hff (cf. 
Eq. (0)), fujj{q) = 2S (j{q = 0) - J{q)Y J{q) is the Fourier transform of the 



exchange integral J^ . Usually the spin wave energies will be smaller by about two 
orders of magnitude than other typical energies of the system as the exchange constant 
J or the Bloch bandwidth W . As a general result the spectral density Ski{E) consists 
of two structures corresponding to special elementary excitation processes of the J, 
electron. There is a rather broad structure built up by "scattering states" which 
result from magnon emission by the original | electron. Thereby the excited electron 
reverses its spin becoming a | electron. Such a process is possible only if there are 
I band states within reach for the original J, electron to land after the spinflip. The 
scattering states therefore occupy the same energy region as the j-DOS (j^). 

There is another possibility for the J. electron to flip its spin. It can also be done by 
a repeated emission and reabsorption of a magnon by the conduction electron resulting 
in a "dressed" particle propagating through the lattice accompanied by a virtual cloud 
of magnons. For not too small positive (negative) J the energy of this "dressed" 
particle lies above (below) the scattering spectrum giving even rise to a bound state, 
i.e. to a quasiparticle with infinite lifetime which we call the "magnetic polaron" . 
Outside the scattering region the polaron peak manifests itself as a (5-function. As soon 
as the peak dips into the scattering part the polaron gets a finite lifetime after which 
it decays into a t electron plus magnon. Figure |l| shows the down-spin quasiparticle 
bandstructure as derived from the respective spectral density as a density plot. The 
degree of blackening is a measure of the spectral density magnitude. Sharp dark lines 
refer to pronounced peaks in the spectral density representing quasiparticles with 
long life-time. For weak coupling \ JS\ < 0.1 scattering processes smear out a little 
bit the "free" dispersion but do not lead to strong deformations. However, already for 
moderate couplings \JS\ > 0.2 one recognizes for some k vectors the appearance of a 
sharp polaron dispersion. For J > (right column) the magnetic polaron is stable on 
the high-energy side of the | spectrum, for J < on the low-energy side. In addition 
the scattering spectrum is clearly visible taking away a great part of the total spectral 
weight. In the antiferromagnetic KLM the magnetic polaron represents the ground 
state configuration For still rather moderate couplings of \JS\ > 0.3 the polaron 
dispersion has split off over the full Brillouin zone. The magnetic polaron has an 
infinite lifetime. It is surprising that even the broad scattering structure is obviously 
bunched together as if it were a rather stable quasiparticle. It is noteworthy to repeat 
that the results of Fig. |l| are exact and free of any uncontrollable approximation. So we 
have to expect these quasiparticle effects in real systems, too. This holds also for the 
quasiparticle density of states (QDOS) plotted in Fig. ^for several exchange couplings 
I J|. According to Eq. (^) the j-QDOS is only rigidly shifted to higher (lower) energies 
for J < (J > 0). Correlation effects appear exclusively in the J. spectrum. For 
\JS\ > 0.25 a band splitting sets in. One of the subbands occupies the same energy 
region as p-\{E) being therefore built up by the mentioned scattering states. In the 
ferromagnetic (J > 0) KLM it is the low energy part of the spectrum containing a 
considerable amount of |-spectral weight. This is not a specialty of the S = ^ case 
or of weak and moderate couplings exhibited in Fig. ^ but holds equivalently, e.g., 
for 5 = 1 and in the strong couphng region {\JS\ > 1.0, see last row in Fig. ||). 
After the band splitting has set in the weights of the two spin-| subbands are close to 
the zero-bandwidth values 2s+i ^'^^ lower, and 23+1 ^'-"^ ^^'^ upper part (J > 0), 
i.e. independent of J. The very often used assumption, when the KLM is applied 




Figure 2. Spin-f (].) density of states of the ferromagnetically saturated 
semiconducter as function of energy as solid (dotted) line. The upper three 
figures of each column show the DOS for a 5 = | system with interaction 
strengths |J| £ {0.1; 0.5; 1.0} corresponding to the spectral densities of Fig. hi. 
The respective lowest figure shows a 5 = | system with |J| = 1.0. The left 
(right) column corresponds to J < (J > 0). 



to the manganites, that the eg electron orients its spin paraUcl to the t2g-S — | 
spin p5| , |2^ , appears with respect to the exact results in Fig. ^ rather questionable. 
In the antiferromagnetic ( J < 0) KLM the low energy quasiparticle subband consists 
of stable polaron states which build the high energy part in the ferromagnetic (J > 0) 
KLM. Although the J > 0-QDOS and the J < 0-QDOS appear as mirror-images the 
implied physics turns out to be rather different. 



3. Dynamical mean-field theory 

3.1. Mapping onto an impurity model 

It is not possible to directly apply the methods of the dynamical mean- field theory [pC| 
to the KLM. One exception is the case of the classical-spin limit (5 oo) 25, 281 , 
thus removing the quantum nature of the spins. The effect of quantum mechanics 
is strongest for S* = -i, but also for S" = |, they dramatically influence the 
spectrum |lj] . The assumption of classical spins for S = | , as done in Refs. |2^ , 
needs therefore careful analysis of the neglected effects. 

Another possibility to derive a dynamical mean-field theory for the Kondo lattice 
model is the fermionization of the localized spin operators as suggested in . This 
approach is, however limited to = ^ systems (see also |5^, Q). The first part 
of our approach will closely follow Ref. [Q, but as discussed below, will then differ 
from the cited reference. 
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The locaUzed spins in the Hamiltonian can be expressed in terms of auxihary 
fermion operators f.^ (ft) i^- 



Here, Tq-o-' represents the Pauh matrices. This is the same transformation that led 
to (Q) where it was applied to the conduction electron spin cr^. The "fermionization" 
according to Eq. (^) implies the introduction of the constraint: 

a a 

The inclusion of this constraint via a Lagrange formalism corresponds to the addition 
of a Hubbard-like interaction term for the /-fermions to the Hamiltonian with the 
interaction constant U'^^^ — > oo The one-particle energy of the /-fermions is 

located at The "fermionizcd KLM" takes the form: 

= ^ .(fc)cLc.. - ^ E s'/' - ^ E flf^^ + ^ E -^-^^ (11) 

k.tj i 1(7 ia 

This Hamiltonian describes a system of two different kinds of = ^ fermions coupled 
by a local spin-exchange interaction. The /-fermions are additionally correlated via 
the Hubbard-type of interaction to prevent double-occupancy. This Hamiltonian 
resembles the periodic Anderson model (PAM) |^ in some way. In fact the only 
difference in the operator structure is the "coupling" between conduction band 
and / state. In the PAM, this is simply a kinetic-energy term ("hybridization"), 
whereas here, the two subsystems are coupled by the spin exchange. Essentially, 
the "fermionizcd KLM" is a rudimentary version of the general multi-band Hubbard 
model with (local) inter-band interaction. For this model, the DMFT is discussed 
by [M- It is straightforward to apply the standard methods of the DMFT to 



map model (11) onto an appropriate impurity model with the corresponding self- 
consistency condition to determine the parameters of the impurity model. This leads 
to the following Hamiltonian for the impurity model (single-site Kondo model, SSKM) 
with fermionizcd /-spins: 

H^-J = E - • s^^^ - ^ E /^/^ + ^ E -'J'-'^l (12) 

with only one "impurity" site for the /-fermions. The spin-exchange interaction acts 
only at this single site denoted by the site-index 0. It is advisable to express the 



conduction band part of Hamiltonian (12) in local space and single out the operators 
referring to the impurity site 0: 

E ^(^)cLcfc<, = J2 '^ij^aCja = (13) 

= ^ToocJ^coa + E^« {4)aCta + cLcoa) + E ^y'^Lcja 

For better readability, we will denote the construction operators at the site with the 
symbol d^r (d^), the on-site energy will be denoted Tqq ~> e^. Finally we introduce a 
unitary transformation which diagonalizes the last term of the second line in Eq. (p^. 
The transformed fermion operators will be denoted aua (^fco-)' hopping T^i Vkd 
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and finally fj{k) for i ^ j. In the context of the DMFT, one does not need 

to know the explicit structure of this unitary transformation since f]{k) and V^d need 
not be known explicitly either. The parameters of the single-site Kondo model will be 
determined by the self-consistency condition of the DMFT (see below). 

A direct solution of model (12) is, to our knowledge not possible. However, it 
can be further simplified: The Hubbard term originating from the constraint (|o|) 
can be eliminated by simply reversing the "fermionization" procedure which led from 
Hamiltonian (^ to (|ri|) . The auxiliary fermion operators get replaced by a local spin- 
i operator at the impurity site, and the constraint (^^ can be dropped. This leads 
to the final version of the single-site Kondo model: 

^^SSKM = Yl ^(^)aLafc'^ + erf 51 ^^'^'^ + Yl ^'^'^ (4afe^ + '^la'^'y) - 'l^^'^'^^ ' S (14) 

k,i7 <y k,a 

This last step, which could be called "de-fermionization" , distinguishes our approach 
from the one used in Ref. and others. This "de-fermionization" ensures the exact 
fulfillment of the constraint (p^), which could for example only be kept on average 

= 1 instead of = 1) in H- 

The parameters determining the conduction band in the Hamiltonian (p^, 
namely ri{k) and Vkd have to be specified according to the DMFT self-consistency 
condition: The local conduction band Green function (cf. (^) should be equivalent to 
the d-operator Green function of the single-site model, ((do-; dj.)) = g'^^\e): 



^GkAE) = Gi^HE) = -— (15) 

where the right equation follows from the formal solution of the equation of motion 
of G^a \E). So instead of the usual definition for the hybridization function, A(i?) ~ 
Sfc E-ij(k) ' ii" i^^s to be determined so that Eq. (|l5|) holds. One will see below that 
the knowledge of A^{E), which can become spin-dependent through this procedure, 
is sufficient to solve the single-site Kondo model (14). Its dispersion fi{k) and the 
hybridization parameter Vkd need not to be determined explicitly. 

The DMFT, i.e. the mapping of the KLM onto the single-site model is, except 
for the limit of infinite spatial dimensions, an approximation, equivalent to the local 
appr oxim ation. In the exactly solvable case of the ferromagnetic semiconductor (cf. 
Sec. 2.3) this is equivalent to neglecting the magnon energies which are assumed to 
be at least one order of magnitude smaller than the energy scales under consideration 
here, e.g. bandwidth or J. 

Next, we will introduce an approximative method to solve the single-site Kondo 
model defined by Hamiltonian (|l3) for an arbitrary hybridization function Aa-{E). 



3.2. Hybridization approximation 

In the following section, we derive an equation of motion-based method to solve the 
single-site Kondo model (|lj). Since in this section we deal only with this model, 
we suppress all subscripts distinguishing between quantities in the lattice and the 
single-site model. 

Starting point is the equation of motion for the d-Grecn function: 

EGi^\E) = h + edGi'\E) + ^ F^dK.; 4)) - ^ E (^-^-i^) + ^-(^)) (16) 

k a 
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where the higher Green functions I^{E) = ((d„S'^;4)) and F„{E) = {{d_^S-'';d^^)), 
corresponding to the Green functions Iim.ja{E) and Fi^jaiE) for the lattice case, are 
introduced. The "mixed" Green function {{uka', rf^)) can be ehminated by investigating 
its equation of motion: 

=T. E^%kf ^'^^^^ ^ ME)Gi^HE) (17) 

thereby defining the hybridization function A{E) = E-i){k) • This yields the final 
equation of motion: 

{E (e, + A(i?))) G^^He) = h-^J2 (^-^'^(^) + ^-(^)) (18) 

Eq. (|r^ ) looks, except for the A(_E) in the prefactor of G^/^E), like the equation of 
motion of the zero-bandwidth limit. The hybridization function A{E) is due to the 
Vfed-term in the Hamiltonian ([l^). This term prohibits an exact solution. In fact, the 
"hybridization" via Vkd is nothing else than the inter-site hopping which was neglected 
in the zero-bandwidth limit. This term will force us to make certain approximations 
in the determination of the higher Green functions on the right-hand side of Eq. (|l8|) . 
We will exemplify this using the /o- (i?)-function. Its equation of motion reads 

{E - ed)h{E) = h{S'') + Vkd{{akaS'-, 4)) - (19) 

fe 

{z„{{d„S'S'-d)j + {{d^^S-'^S^'-dl)) - z„{{d^dU-aS-^;di))) 

where on the one hand, higher "impurity-site" Green functions are introduced, but on 
the other hand also a higher "mixed" Green function, {{akaS^; dj.)). In analogy to the 
one-particle mixed Green function, Eq. (|l7|), we use the following substitution: 

J2VM{{akaS^;4)) ^ AiE){{d,S^;dl)) ^ A{E)I,iE) (20) 

k 

The justification of this procedure can be found in analogy to the "self-energy 
substitution" known from other approximation methods for the KLM (see e.g. ]T^) 
by inspecting the spectral representations of the relevant Green functions. This reveals 
that all of them have the same single-particle poles, they differ only in the respective 
weights given by matrix elements of the type {En\A\Em) ■ Here |i?„) are the energy 
eigenstates and A represents one of the relevant operators: da, dka, d^S^ or a^aS^ . 
From Eq. ( p7[ ) follows that the hybridization function accounts for the differences of 
the matrix elements of the first two operators {d^ and aka)- Assuming that these 
differences are almost equal to those of the matrix elements built up by the latter two 
operators leads to the hybridization approximation (^). 

This is the only approximation necessary to decouple the hierarchy of equations 
of motion. In the case of 5* = ^ spins, there are only 6 different "impurity-site" 
Green functions, whose equations of motion form a closed set of equations. Besides 
the already introduced Green functions G^'a^ (E) , 1^ (E) and F^ (E) , these are 

Fi'\E) = {{d^adl,daS^-;dl)) 

F^'\E) = {{d^adUaS-'';di)) (21) 
Da{E) ={(d^adlada;di}) 
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Several expectation values, introduced into the theory via the inhomogeneities of the 
equations of motion can be expressed in terms of the above-mentioned Green functions, 
a self-consistent solution has to be found. 

At this point, let us comment on the reliability of this approximation. Although 
only one approximation enters our decoupling procedure (cf. Eq. (pO|)), it still has 
to be seen as an uncontrollable approximation meaning that there is no true small 
parameter. There are two non-trivial limiting cases where the replacement ( pO[ ) 
becom es ex act: The first is the limit Vkd 0, representing the zero-bandwidth KLM 
of Sec. |2.2| . But already a small, but finite bandwidth in the KLM could lead to any 
(unknown) expression for A(E) implying that we can not necessarily assume Vkd to 
be small any more. The second limit is the "classical spin" limit where the "spin 
variable" 5' has no operator properties any more. Here, the replacement of Eq. ( po|) 
reduces to (0). However since we are interested in the general case with finite S 
and bandwidth, we have to confirm the trustworthiness of this approximation by a 
comparison with other well-tested methods. 



4. Other approximation methods 

We are now going to introduce three further approximation methods for the KLM. 
These are known from literature, their strengths and weaknesses have been identified. 
The three methods differ substantially with respect to the theoretical assumptions 
made for an approximate solution of the KLM. Common features following from these 
procedures and the above-introduced DMFT scheme should then give some credit of 
reliability , in particular when they additionally fit the exact limiting cases discussed 
in Sec. |. 

4-1- Second- order perturbation theory 

An application of the usual diagrammatic perturbation theory to the Kondo model 
cannot be performed due to the absence of Wick's theorem. Only for low temperatures 
(spin- wave approximation) or in the classical-spin limit |5^, this method is 
applicable. The projection operator formalism of Mori |4^, ^ is better suited for 
the Kondo model. It has been used successfully to describe correlation effects in 
the Hubbard model in the weak coupling regime JSTf . The general formula for the 
second-order contribution T,^^^{E) can be found there (Eq. (3.12) in Ref. |5^). To 
allow for a better comparison with the other approximations in this paper, we further 
approximate the self-energy taking only k averaged occupation numbers into account 
(local approximation) . 



4-2. Self- consistent CPA 

Next, we want to introduce a modification of the well-known "coherent potential 
approximation" (CPA) for the KLM. The CPA is a standard many-body approach 
that starts from a fictitious alloy in analogy to the interacting particle system. Starting 
point may be the zero bandwidth limit of Sec. 2.2. We think of a four-component alloy 
each constituent of which is characterized by one of the energy levels Ei in Eq. (||). 
The spectral weights (^) are then to be interpreted as the "concentrations" of 
the alloy components as seen by a propagating cr-electron. The fictitious alloy for a 
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cr-electron is built up by the local moments and by the frozen (— cr) electrons. The 
CPA-selfenergy of the cr-electron is found by the well-known formula [5^ . 

As a consequence of the single-site aspect of the CPA the resulting selfenergy is wave- 
vector independent. According to Eq. the "concentrations" a^^ depend on a sum 
of the "higher" correlation functions lu^iaiE) and Fa^icriE), which can rigorously be 
expressed by the single-electron Green function. 

— 1 1 f^°° 

7, + z,A, dEf^{E){E-e{k))lmGka{E) (23) 

The shortcomings of the CPA-procedure lie on hand. The one is the same as that 
in the conventional alloy analogy of the Hubbard model, namely the assumption of 
frozen (— cr) electrons. This is partially removed by our proposed modification of the 
standard CPA procedure for the KLM, namely the sclfconsistent calculation of the 
higher correlation functions 7o- and A^- via Eq. (|2^) as well as the band-occupation 
via the spectral theorem for the one-electron Green function. Maybe even more 
serious in the case of the KLM is the blocking of repeated spin exchange with the 
local moment system. Magnon emission or absorption is not involved. So we cannot 
expect that the CPA-treatment correctly reproduces the exact limiting case of Sec. 2.3. 
However, some general information about the quasiparticle bandstructure might be 
possible, in particular in the strong coupling ("split band") regime. By construction 
the method yields the correct zero-bandwidth limit. 



4-3. The moment- conserving decoupling approach (MOD A) 

A Green function method which takes the spin dynamics correctly into account has 
been proposed in Ref. iQ. For details about this approach, we refer the read to the 
cited paper, here we summarize the result shortly 

This decoupling approach yields finally a selfenergy of the following structure: 

^ka{E) = -]^JzAS'-) + ^J^Dua{E) (24) 

The first term is linear in the coupling J and proportional to the 4/ magnetization 
{S^). It just represents the result of a mean-field approximation being correct in the 
weak-coupling limit. The second term in ( p^ ) contains all the spin exchange processes 
which may happen. It is a complicated functional of the selfenergy itself. So (|2^) is 
not at all an analytical solution but an implicit equation for the selfenergy. We do 
not present here the lengthy expression for Dk„{E) referring the reader for further 
details to Ref. ||lj|. It should be mentioned, however, that Dka{E) contains several 
expectation values which must be fixed to get a self-consistent solution. No problems 
arise with the mixed correlation functions 7^ = {Sfcl_^Cia) and A^ — {S^Uia). They 
can rigorously be expressed by the spin flip function Fn^i„{E) and the Ising function 
Iii,ia{E) defined previously. Both functions arc already involved in the procedure, so 
that no further approximations are necessary to fix and A^.. For pure local-moment 
correlations such as (S^), (SfS^^),. . . , however, a special treatment is necessary, e.g. 
as described in Ref. [|14|. 




Figure 3. Conduction band density of states (DOS) obtained via the DMFT, 
MCDA, CPA and SOFT (from left to right) for different values of J as function of 
energy. The temperature is set to T = 0, therefore (5^> = 0.5 (CPA: (S^) = 0.33, 
see text) and the electron density n = 0.5, the chemical potential is at /i = 0; 
Solid line: spin-|, dotted line spin- J, DOS 



5. Results 

In the following section, we discuss the results obtained for the DMFT and the three 
approximation schemes of Sec. ^for the ferromagnetic Kondo lattice model with S = 
Since we are interested in the reaction of the conduction band due to the magnetic 
order of the spin system, we have not calculated the latter self-consistently. Instead, 
we have simulated the magnetic order by determining (S^) using a Brillouin function. 
Temperatures are given in units of Tc- Within the CPA, our choice of (S^) can lead 
to unphysical results. Namely in the case of {S^) S, some of the weights 
can become negative. There is an upper bound for (5^) [Q. We therefore had to 
limit (5^) to (S^) < 0.33 for some of the CPA calculations. Within the DMFT 
calculations, we experienced severe numerical problems which forced us to introduce 
a further approximation: "mean- field" -decoupling the Fa^\E) and Fi'^\E) functions 
simplifies the system of equations of motion: 

Fi'\E) « (1 - (n_,))r.(i?) + {S^)D4E) + +((1 - n^^)S^)Gi'\E) 

F^^\E) « MF,{E) {S--dU-a)GifHE) (25) 

All DMFT-results presented below were obtained using the hybridization approxi- 
mation in combination with this unrestricted-mean-field decoupling of i^i^' (E) and 

f^^\e). 

In all calculations, the conduction band is described by a tight-binding DOS for a 
simple-cubic lattice structure of unit width {W — leV). The Curie temperature 
is taken as Tc — 250K. 



Quantum effects 



. Kondo lattice model 



15 



DM FT 



MCDA 





CPA 



J=0.1 




SOFT 








J=0.2 




J=0.3 






J=0.4 




Figure 4. Same as in Fig. H, but for T = Tc, therefore (S^) = 0. 



The quasiparticle densities of states (DOS) for quarter-filling and different values 
of J are plotted in Figs. ^ and || for T = and T = Tc, respectively. As indicated, 
the columns correspond to DMFT, MCDA, CPA and SOFT, respectively (from left 
to right). 

We will begin the discussion with the DMFT results shown in the left column 
of Fig. ^. A small value of J leads to a spin-dependent shift of the spin-f and J, 
DOS, as one would obtain by a simple mean-field decoupling, i.e. by replacing S"^ 
by its mean value (5"^) in Hamiltonian (|l|). With increasing J, the DOS show some 
striking correlation effects: first a broadening, later the onset of a splitting of the band 
can be observed. Whereas in general, the correlation effects are stronger for spin J, 
(indicated by a stronger quasiparticle damping), the splitting is, in contrast to the 
other two methods, more pronounced in the spin-f DOS. However, here the split-of 
(upper) peak has much less spectral weight than the original peak which, except for 
a band-narrowing still resembles strongly the free conduction band DOS. 

The picture is very similar in the MCDA. Again, a mean-field like spin-dependent 
bandshift is observed for small J. On increasing J, the spin-| DOS broadens, and a 
two-peak structure emerges. The spin-f DOS remains, except for a small tail at its 
upper edge, unchanged. This behavior can easily be understood by comparing with the 
special case of the ferromagnetically saturated semiconductor as discussion in Sec. 2.3 
since the MCDA develops continously into this special case for n — > and T — > 0. 
The spin- J, DOS splits into a scattering part (low energies) and the polaron-like part 
at higher energies above the conduction band. Unlike the n = case, however, there 
are some, but weak modification in the | DOS, namely the above mentioned tail at its 
upper edge. This can be interpreted as a scattering contribution for spin-t electrons. 
The origin of this is the finite number of spin- J, electrons. 

Also the CPA DOS show a remarkably similar picture as the previously discussed 
theories. Again, there is the mean-field shift for small J. For larger values of J, the 
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splitting of the spin-| DOS sets in similarly to the other two methods. Here, this can 
be contributed to the two single-occupancy quasiparticle energies Ei and E2 known 
from the zero-bandwidth limit (cf. Sec. 2.2). For n 7^ a third quasiparticle energy 
carries non-vanishing weight, namely the E4 peak. This is located between the two 
other peaks since we switched off the conduction band Coulomb interaction U (cf. 
Eq. (§). For J = 0.4, the appearance of this third peak in between the two main 
peaks is vaguely visible. A remarkable difference to the DMFT and MCDA results is 
the fact, that the spin-f and scattering part of the spin- J. DOS do not cover the same 
energy range. This shortcoming of t he CPA is due to the neglection of spin-dynamics 
(magnons) as mentioned in Sec. O (cf. Ref. po|] ). 

The SOFT results are, for the plotted values of J in the weak- and intermediate 
coupling regime, not too far away from the other results. For small J, where the 
SOFT becomes by definition reliable, the mean-field shift as in the other methods is 
clearly observable. Similar to the other methods, the deformation of the spin- J, DOS 
is much stronger than that of the spin-f DOS. However, for larger J, the SOFT never 
shows a true band splitting, its range of validity is certainly restricted. 

For T = Tc, where spin-symmetry is re-established, the DMFT, MCDA and CFA 
give a similar overall picture. The resulting DOS shows, already for J > 0.3, a two- 
peak structure (DMFT and MCDA), in the CPA a third peak is dimly noticeable. 
The most remarkable observation is the near-coincidence of the DMFT and MCDA 
results for all J values. The transition from the (clearly distinct) T = DOS to the 
(resembling) T — Tc DOS is continuous for both theories. The SOPT, however, has 
to be seen as a complete failure for a paramagnetic system already for relatively small 
J « 0.3. This can already be read of the formula for calculating the self-energy, which 
is more or less trivial. The corresponding results cannot be connected to any of the 
exactly solvable limiting cases (cf. Sec. 4.1). 

A sound-standing interpretation of the observations is not difficult since the 
MCDA can be traced back to the exactly solved limiting case of the ferromagnetically 
saturated semiconductor (cf. Sec. |2.3D and the CPA to that of the zero-bandwidth 
limit discussed in Sec. The SOPT, of course, becomes reliable for small J. 

Although all four presented methods are of approximate nature and their results 
do, at least for intermediate-to-large values of J, differ, some common properties 
emerge: For small J, the behavior is genuine mean-field like, a bandshift proportional 
to ±J{S^) is observed. For larger J, the onset of a band splitting occurs. At T = 
this primarily affects the spin-| DOS. This fact is understandable by examination of 
the ferromagnetically saturated semiconductor where the band splitting was discussed 
in terms of a scattering band and the magnetic polaron. In that case the spin-t 
DOS remains, except for a simple shift, unaffected by the interaction simply because 
spin-flip of spin-f electrons is suppressed in this case. Now for finite n, this is only 
approximately true. The spin- J, DOS still shows strong correlation effects, but also 
the spin-| DOS is affected due to the finite number of spin- J, electrons in the system. 
This manifests itself differently in the various methods: in the MCDA, a tail is seen at 
the upper edge of the spin-t DOS, in the CPA, a shoulder develops, and in the DMFT 
approach, a band splitting is observable. From this we conclude that correlation 
effects are much more pronounced in the DMFT than in the other two theories. At 
T = Tc, the dip indicating the onset of the band splitting is still existing for all but 
the SOPT method. In CPA, MCDA and the DMFT results this splitting is of similar 
size, which can easily be read off the CPA where it is simply given by the difference 
of the respective energies from the zero-bandwidth limit (see Eq. (0)). The two most 
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pronounced peaks correspond to the single-occupancy quasiparticle energies Ei and 
E2, the band sphtting is therefore approximately AE = J{S + i). 

This is in contrast to the results obtained by dynamical mean-field theory for the 
S Qo KLM (KLM with classical spins). The emerging picture for classical spins is 
the following At T = 0, the DOS is characterized by a mean-field like "Zeeman" 
splitting between the bands of both spin directions. With increasing temperature, at 
each of the respective spin-| or J, band, spectral weight of the opposite spin direction 
appears, until finally at T = Tc, the system becomes paramagnetic. The splitting into 
two subbands separated by AE = JS stays constant. Comparing our results with the 
5 — * 00 ones, phenomenologically the DOS's show completely different characteristics 
at T = 0, especially in the spin-| band. For T = Tc, the DOS's look quite similar. 
However, the physics behind the scenes turn out to be completely different as can 
be seen, e.g., by the inconsistent size of the band splitting AE. Whereas in the 
classical-spin limit, the splitting is always due to a mean- field like "on-site Zeeman 
splitting" , in the case of quantum spins, the various elementary excitations as discussed 
in the case of the ferromagnetically saturated semiconductor, are responsible for the 
sub-band structure. For the ferromagnetically saturated system, the neglection of 
spin-flip processes and magnons in the S* — > cx) calculation lead to the picture of a 
half- metal |2^. This does not apply for any finite value of S and is therefore clearly 
an artifact of the 5* — > 00 limit. 

Let us shortly remark on the situation with anti-ferromagnetic coupling (J < 
0) |]60|| . All four approximation methods presented here do not show any signs of 
Kondo screening. An investigation of the special low-temperature physics of the model 
with J < is therefore not possible. However, some remarks about the behavior of 
the model for T > Tk can be made: In general the excitation spectra are broader than 
in the J > case. The resulting exchange splitting of AE « J{S + 1) is also larger 
than for J > which can already be seen in the zero-bandwidth limit. This is again 
in sharp contrast to the S* = 00 results, where the size of the splitting is independent 
of the sign of J. 

6. Conclusions 

In this paper, we investigated the Kondo lattice model (KLM) focusing on the model 
with positive (ferromagnetic) exchange constant J. We discussed two exactly solvable, 
but nevertheless non-trivial limiting cases as well as four different approximation 
methods. The results obtained by these methods, maybe except for the perturbation 
theory, compare generally reasonably well. Sometimes even nearly perfect matches 
occur (cf. Fig. ^). In general, the differences between the methods are smaller for the 
paramagnetic than the ferromagnetic system. 

From the comparison of common features of these approximation methods in 
combination with the exact results, the following picture emerges for the quasiparticle 
structure of the ferromagnetic Kondo lattice model: For small | J|, a mean-field like 
shift of t and J, DOS is visible. Already for intermediate coupling strengths, a band 
splitting of size AE = J{S + ^) occurs. The relevant energy scale for the splitting 
is JS « 0.3. In the case of ferromagnetic saturation (T = 0), the splitting is more 
pronounced in the | than in the f DOS. This band-splitting should not be confused 
with the splitting found in the limit of "classical spins" {S — * 00). There, the 
splitting is simply due to a mean-field like Zeeman-splitting and therefore of the size 
AE — JS. Contrary to that, our results clearly show that the two emerging subbands 
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can be traced back to the two elementary excitations known from the hmit of the 
ferromagnetically saturated semiconductor (see above). The magnetic saturation of 
the spin-system suppresses these processes for spin-l electrons which explains the 
stronger footprint of the correlations in the spin- J, DOS. For finite temperatures, the 
magnetic polaron generally remains a well defined quasiparticle, represented by a 
rather sharp Lorentzian peak in the spectral density. Now the spin-| electrons can also 
participate in spin-exchange processes since the localized spins are not fully aligned 
any more. The spin-symmetric DOS at T > Tc also show the characteristic splitting. 
Our results also confirm the fundamental differences between the J > and the J < 
case of the Kondo lattice model: It can be read of both exactly solvable limiting cases 
that the ground state will be different depending on the sign of J. The approximative 
approaches of Sees. || and ^ do not allow a deeper investigation of the model with 
J < due to the inability to reproduce the special low-temperature properties of that 
model ("Kondo physics"). Another important conclusion from our calculations can be 
drawn: There is always finite spin- J, spectral weight in the region of the spin-f DOS. 
This spectral weight does not disappear in the limit J ^ oo but only in the limit 
S* — > oo ("classical spins"). For the ferromagnetic KLM (J > 0), this implies that the 
KLM for 5* = I and large J, corresponding to the situation found in manganites, will 
not be a half-metal for T — 0, contrary to the predictions of 5" — > cx) calculations. 
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